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ABSTRACT 

We describe the effect that new atomic calculations, including fuUy-relativistic R-matrix calculations 
of collisional excitation rates and level-specific dielectronic and radiative recombination rates, have 
on line ratios from the astrophysically significant ion Ne IX. The new excitation rates systematically 
change some predicted Ne IX line ratios by 25% at temperatures at or below the temperature of 
maximum emissivity (4 x 10^ K), while the new recombination rates lead to systematic changes at 
higher temperatures. The new line ratios are shown to agree with observations of Capella and ct^ CrB 
significantly better than older line ratios, showing that 25-30% accuracy in atomic rates is inadequate 
for high-resolution X-ray observations from existing spectrometers. 
Subject headings: atomic data, line: formation. X-rays: stars 
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1. INTRODUCTION 

X-ray spectral models for collisionally ionized plas- 
mas have improved steadily over the past 30 years. 
As spectral lines dominate the emission from an X- 
ray emitting colli sional plasm a over a wide range of 
temperatures (Co x fc Tucker I IT969). a relatively large 
number of lines must be included to create an ac- 
curate model. However, most early X-ray observa- 
tions had limited spectral resolution, R = E/AE ~ 
1 — 20, and approximate calcula tions of individual 
transitions w ere ge n erally adequate (iRavmon d & Smi^ 
119771 iMewd 119721 : iLandini fc M onsignori Fossil 11972^ 
As laboratory and observational methods have im- 
proved, more compl et e calculations have been produced 
(iSmith et al.l 120011: iKaastra. Mewe fc NieuwenhuiizenI 
119961: iDere et al.l ll997D. although questions remain as 
to whether or not these are sufficiently accurate for as- 
trophysical applications. With the launches of Chandra 
and XMM-Newton, individual emission lines from a wide 
range of elements are regularly resolved and used in diag- 
nostics, e specially from the grat i ng instruments on thes e 
missions (iCanizares et al.l 120001: iBrinkman et al.ll2001[ ). 
The strongest lines are from He-like and H-like ions and 
the Fe L-shell ions. Since both terrestrial and astrophys- 
ical plasmas are dominated over a broad temperature 
range by emission from the complex Fe L-shell ions, a 
collabor ation of atomic phys icists established the Iron 
Project (jHummer et al.|[l993D to perform high-accuracy 
R-matrix calculations of these ions. However, no equiva- 
lent systematic calculation of line emissivities for all He- 
like or H-like ions exists. We present here an initial study 
of Ne IX comparing new high accuracy calculations of 
collisional and radiative transition rates to a number of 
existing results. 

We began this project with He-like ions because they 
provide powerful diagnostics on the physical conditions 
of the plasma. Neon is an obvious choice as it is an 
astrophysically abundant element not affected by dust 
depletion, with emission lines around 1 keV that can be 
easily detected by both Chandra and XMM-Newton. Fur- 
thermore, experimental measurements of Ne IX spectra 
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Fig. 1. — Energy level diagram for the n = 1,2 levels of Ne IX, 
with radiative transitions marked in colors that show the radiative 
transitio n rate (in s~ ^). The strong w,x,y, and z transitions are 
marked IIGabriellll972l ). The 2- photon transition is shown with two 
extra arrows. 



are available to test the models () WargelinI 1 1 9931 ) . In col- 
lisional equilibrium models, at least 1% of neon is He- 
like in plasmas with temper atures between 5x10^ and 
10^ K (jMazzotta et al.lll998t) . and thus Ne IX emission is 
present in numerous types of astrophysical sources. 

The transitions from the n=2 He-like energy levels 
are particularly favorable for use as diagnostics, because 
their rates vary widely. Figure [1] shows the fine-structure 
levels arising from n=2 together with the transitions: w, 
X, y, and z, with eight orders of magnitude variation be- 
tween the w (resonance line) and the z (strongly forbid- 
den line) transition rates. The emission lines from these 
transitions in He-like ions have been used in a number of 
ways, including not only in the well-known G-ratio (G = 
{x + y + z) /w) for temperature and i onization state and 
R-rat io {R = z/{x + y)) for density (jGabriel fc Jordanl 
Il969f ). describe d below, but also as d iagnostics for 
photoionization (Porquet & Dubau '2000^ , photoexcita- 
tion (,Mewc fc Schriiver, ,1978. : .Kalin et al.. ,2001.) . and 
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non-M axwellian velocity distributions ()Gabriel et alj 
Il991f ). Early observations of He-like ion lines, in- 
cluding Ne IX, in solar active region and flare X- 
ray spectra have provided general diagnostics for den- 
sity and temperature, but with complications due to 
line blending, resonance scatter ing, and atomic data 
poschek 1 119901: iMcKenzid I1987D . More recently, for 
stellar coronae, Ne IX emission lin es have been used 
to measure relative Ne abundances (iDrake et al.l 12^0011 : 
iHucncmocrder. Canizares fc SchuI3 120011 ) and optical 
depths (Ne ss ct al . 2003a) , as well as temperatures and 
densit ies (jAvres et alJ200ll:lOsten et al.ll2003l : lNess et alJ 
l2003bD . Measurements of changing line ratios dur- 
ing stellar flares are used to estimat e the temperature 
and density in the hea ted plasma (Giidel et alJ 120021 : 
Ivan den Besselaar et al. 2003; Raasscn ct al. 2003). 

He-like emission line ratios have been calcula ted by 
a num ber of authors, including ,Bautista & : Kallman I 
(|2000[ ). who considered both co llisional and photoion- 
ize d plasmas. iSmith et al.l (|2001[ ). who focused on O VII, 
by iPorguet et al.l ( 2001|) . computing mod els for C V, 
N VI O VII, Ne IX, Mg XI, and Si XIII. ISmith et all 
()2001h showed that including higher-n states, espe- 
cially level-resolved dielectronic recombination to high- 
n states, could significantly aff ect the G-ratio at higher 
energies. IPorguet et al.l ()200lD included the effect of a 
radiation field on the forbidden and intercombination 
lines, but these calculations were based on distorted 
wave (DW) electro n -impa ct excitation rate calculations. 
iZhang fc Sampsoiil ()1987n considered resonance effects 
manually added to a DW calculati on, but only for th e 
n = 2 states. These were included bv lSmith et al.l (|200in . 
with little resulting effect. 

Despite advances in atomic calculation meth- 
ods and the frequent use of He-like diagno stics, 
Bautista & KallmaT] (|2000l ). ISmith et all (|200lD and 
Porguet et al.. (■20011 ) all used rates taken from the DW 



calcul ations, primarily from iSampson. Goett &: Clark! 
(|1983f ). Although some modern R-mat rix calculation s 
are now available ([Bautis ta,, 2003; Delaha ve et al.l[2006l) . 
the compilation by S ampson. Goett & Clark! ( 19831 ) is 
one of the few complete collections of He-like ion rates 
that includes all ions of astrophysical interest at X-ray 
energies. A signal that these data may not be adequate, 
however, came from a detailed study of Ne IX lines 
from Capella, where iNess et all (|2003b( ) found that the 
G-ratio diagnostic suggested a temperature less than 
2 X 10^ K, despite the well-known and strong peak in 
Capella's emission measure distrib ution at 6 x lOfK. 
A similar discrepancy was found by iTesta et "all (|2004D 
from the G -ratios of O VII and Mg XI. Meanwhile, 
lOsten et al.l ([20 03) inferred a low electron density for 
(T^ CrB from the Ne IX R-ratio diagnostic, despite find- 
ing larger values from other diagnostics. These problems 
suggested to us that resonances, which tend to increase 
the collisional excitation rate of non-dipole transitions at 
low energies, could be significant. We knew these were 
not included in t he existing line ratio calculations of 
ISmith et al.[ (|200H) that used the DW method. However, 
the R-matrix method will automatically generate all 
relevant resonances, and can be extended by carrying 
out a DW calculation for larger principal quantum n 
as well. A number of such calculations exist, including 



iBautista I ([2003[ ). lDelahave et al.l (|2006D . and lChen et"all 
(|2006D . For this work, we used the R-matrix calculation 
of Chen et al. (2006) for Ne IX, which include the 
largest number of states (up to n = 5). 

We briefly describe the astrophysical model in ij2l in- 
cluding the R-matrix calculation and our sources for 
other data used in the spectral calculation. We com- 
pare our results to previous work in Sj3l and discuss the 
effects on common spectroscopic diagnostics in []4| 

2. METHOD 

Once all the necessary atomic data have been com- 
piled, calculating line strengths for collisional plasmas is 
straightforward with the APEC code (S mith et al. 2001 ). 
APEC uses an explicit rate matrix formulation includ- 
ing collisional and radiative transitions between levels, 
as well as dielectronic and radiative recombination from 
the next higher ionization state to the ground and ex- 
cited levels of the ion under consideration. The APEC 
code does not contain atomic data within the code itself, 
but instead reads from a database of FITS files that can 
be modified as needed. This makes comparisons between 
different atomic calculations easy, as any changes in the 
database are entirely independent of the code. 

2.1. Collisional Excitation Rate Calculation 

For n < 5 , we u se the collisional rate coefficients 
of iChen et al.l (|2006[ ) . These electron collisional exci- 
tation calculations use a f ully- relativis t ic clo se-coupling 
approach. As discussed in iChen et al.l (|2006[ ). the addi- 
tional resonance excitation included in this calculation 
brings the model G-ratio values into much better agree- 
ment with the laboratory results measured at the LLNL 
Electron-Beam Ion Trap (EBIT) experiment ([Wargelinl 
[1991. 

ISmith et al.l (|2001[ ) showed that high-ri transitions (up 
to n = 10) can be significant to the high-temperature val- 
ues of the G-ratio. The variation appears to be largely 
due to dielectronic recombination (see H2.S^ into high-n 
states, but for completeness we added collisional excita- 
tion rates for levels 6 < n < 10 calculated with the DW 
FAC code (IGull2003D . 

2.2. Radiative recombination 

APEC initially used photoionization rates to calculate 
the rate of radiative recombination to ground and excited 
states of Ne IX, via the principle of det ailed balance, pri- 
marily from lVerner fc YakovlevI (119951) for ground-state 
photoionization and iClark. Cowa n fc Bobrowicz ' (1 19861 ) 
for excited-state results. Recently Badncll (2006aj) cal- 
culated level-separated radiative recombination rates for 
all atoms with Z < 30 and ion stages up to and including 
Na-like ions recombining to Mg-like ions. These data are 
significantly improved from the original rates in APED, 
and they extend to higher individual n l evels (n < 8 com- 
pared to n < 5). We therefore use the iBadnelf (|2006aD 
calculation in this paper. We n ote that a com parison 
of the existing APEC results and iBadnelll (|2006al ) shows 
that the two calculations agree to better than 30% for 
plasmas with temperatures between lO'' and 10^ K for 
the rates to the individual n = 2 levels and the total 
recombination rate to Ne IX. 
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2.3. Dielectronic Recombination (DR) 

Dielectronic satellite recombination involves two elec- 
tron transitions. The incoming free electron, with energy 
Ei , recombines into an excited level with energy Ej of ion 
/ while also exciting a bound electron into energy level 
Ek, leaving the ion in a doubly-excited state. Unlike 
radiative recombination, this is a resonant process since 
Ei = Ej + E]^ . The doubly-excited ion then relaxes when 
one electron transitions to a lower state (often a state 
with the lowest unoccupied orbital) and either ejects the 
second electron or emits a "satellite line" photon. This 
latter case is the DR process, and the photon emitted 
appears as a red-shifted "satellite" of the transition that 
would occur without the presence of the second excited 
electron. The ion is then left in an excited state, which 
usually decays via a radiative transition. 

All DR calculations include the explicit generation of 
many different satellite transitions; depending upon the 
configurations included, the number of transitions can 
range into the millions or more. Each transition has a 
temperature-dependent rate that generates a particular 
DR satellite line and leaves the ion in a specific excited 
state. In practice the DR rates are summed either to 
generate the total DR excitation rate for each excited 
state (eliminating information about the satellite lines), 
or to obtain the total DR rate for the ion (eliminating 
information about DR excitation). Of course correctly 
calculating the output spectrum requires individual tran- 
sition data, which are difficult to obtain. At sufficiently 
high resolution, the satellite lines can be separated from 
the primary transitions, so for this project we ignore the 
satellite lines and use DR data that have been summed 
to create level-specific recombination rates. However, as 
noted bv iBautista fc Kallman I (|2000( ). the G-ratio mea- 
sured with spectrometers with low- or medium-resolution 
may be affected by blending of the primary lines with DR 
satellites. For the DR data, we use rates from the recent 
calculation of Badnell (2006bJ3 , which include transi- 
tions up to 71 = 8 for recombinations to Ne IX. 

3. RESULTS 

We are now able to compare line ratio results from a 
number of different calculations. We used four different 
calculations for this paper. The Full model uses the 
complete data set, as described in fJJl while the Full (No 
Rec) model uses the complete data set excluding radia- 
tive and dielectronic recombination. The ATOMDB 
vl.3.1 model uses the original ATOMDB vl.3.1 
database (iSmith et al.ll200lD. This model uses e lectro n 
excitation rates from iSampson. Goett fc ClarkI (|1983l ). 
radiative recombination ra tes derived from the partial 
photq ionization rates of iClark. Cowan &: Bobrowicz I 
(119861). dielectronic recom bination rates from 
iVainshtein fc Safronoval ()1978D and iSafronova et al.l 
and radiative transition rates from a number 
iese et all 11996': 'Derevianko & Johnson ' 
T977r 
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Fig. 2. — Ne IX G-ratio as a function of temperature for four dif- 
ferent calculations, at a density of 1 cm""^. The starred points show 
data from EBIT measurements folded through a Maxwellian distri- 
bution. The EBIT contained no Ne+9 during the experiment, and 
so the EBIT points (including an estimated 8% error) should be 
comp ared to the red curve as no recombination occurred llWareelinI 
[T99a i: see lChen etan lf200^ ') for more details. 

FAC model uses atomi c data tak en exclusively from the 
Flexible Atomic Code (|Gul 120031 ) and does not include 
recombination. These four cases allow discrimination 
amongst changes due to the R-matrix calculation alone, 
those due to a larger transition matrix, and those due 
to increased consideration of recombination effects. We 
postpone to a later paper consideration of effects due 
to innershell ionization from Ne+^ due to the small 
range of temperatures with both significant Ne"*"^ ions 
and detectable Ne IX lines. In all cases we assume 
collisional ionization equilibrium, using the tables of 
'Mazz otta etall (jl998D . 

We focus on two line ratios commonly used for He-like 
ions, the already-described G-ratio, which is a tempera- 
ture and ionization state diagnostic, and the R-ratio, a 
density diagnostic. In Figure [H we show the G-ratio as 
a function of temperature for our four different models. 
A number of effects are immediately apparent: (1) The 
R-niatrix results create a systematically larger G-ratio, 
(2) the FAC results (which do not include resonances, 
but do extend to n < 10) show a more complex behav- 
ior which is generally closer to the current ATOMDB 
results, and (3) including recombination systematically 
changes the result for temperatures above ~ 3 x 10^ K. 
The relative difference between the Full results and those 
in ATOMDB vl.3.1 is an almost constant 25%. 

Figure [3] shows the R-ratio as a function of electron 
density at two different temperatures and as a function 
of temperature (at an electron density of 1 cm~^). As the 
density increases, the slow radiative decay of the ls2s^So 
to the ground state (the z line) is eventually matched by 
electron excitation to the ls2p'^Pi^2 levels, which decay 
to the ground state (creating the x and y lines). For 
a given temperature, the models are all quite similar at 
high densities, but diverge from each other at lower den- 
sities, reflecting in part the additional resonance contri- 
bution in the Full and Full (No Rec) models. Also 
at the lower densities, the Full and Full (No Rec) are 

similar for the lower temperature case (where very lit- 

www-cfa dc.phy.ornl.gov/data-and-cod cs/aur ost/aurost-recomb/ hpii^e7^]:)^a neon is present) but differ significantly at the 
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higher temperature due to recombination. Above the 
low-density hmit, as the R-ratio begins to dechne rapidly, 
the different models lead to a factor of two difference in 
the inferred density. Figure |3[Right] shows that recom- 
bination produces a significant increase in the R-ratio as 
the temperature increases. However, at lower tempera- 
tures the biggest differences among the models are due to 
differences between the R-matrix and DW calculations. 



4. DISCUSSION AND CONCLUSIONS 

The original goal of the project was to determine 
if the unexp ected results found by observers such as 
lOsten et al.r()2003D and INess etHI (|2003bf ) were due to 
astrophysics or simply due to inadequate atomic models. 
In both cases, the authors noted the issue but were un- 
able to distinguish between the two possibilities, and un- 
fortunately, our results confirm that there is no simple or 
quick method to select between them. Our results show 
a difference between the DW and R-matrix of ~ 25%, 
within the oft-quoted error of "approximately 30%" for 
atomic calculations. They also demonstrate that this 
30% inaccuracy is inadequate for use with bright lines 
observed by current X-ray satellites. 

In Figure |4[Left] we show the Ne IX G-ratio in the 
same format as Figure 11 of INess et all mom . The 
horizontal dashed line shows the measured value from 
Cap ella (with errors shown as dotted lines). As noted 
bv Ness et al.l (|2003bf ). the result from ATOMDB vl.3.1 
barely agrees with the observed ratio, and only at a very 
low temperature. Our new calculation predicts a tem- 
perature in the range logT ^ 6.3 — 6.6. Differential 
Emission Measure (DEM) fits of Capella using X-ray 
and EUV spectr a have consistently show n a strong peak 
at logT w 6.8 (|Brickhouse et al.l [19951) . although with 
some dispersion and emission at other temperatures. The 
agreement between our results and laboratory measure- 
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ments (|Chen et al.ll2006[ ) suggests that the remaining dif- 
ference between the G-ratio and DEM temperatures may 
well be due to the plasma conditions, and not to errors 
in atomic physics. 

Although the G-ratio is useful as a temperature di- 
agnostic, temperature can also be measured from the 
bremsstrahlung continuum or from broad spectral fits. 
However, density diagnostics require specific line ratios, 
and the strong lines emitted by He-like ions make the 
R-ratio both useful and usable. Figure |4[Right] shows 
the ratio for all four calculations, done at two tempera- 
tures, logT = 6.4 and 6.8. These were chosen since they 
cover the range from the best-fit temperature using the 
G-ratio and from DEM fits for both Capella ([Ncss et al] 
2003bD. ari d, coincidentally, cr^ CrB (lOsten et al.ii2003l ). 
Ness et all (2G03b) noted that their value for the R-ratio 
in Capella (3.97 ± 0.67, shown as the yellow box) is only 
marginally in agreement with the ATOMDB vl.3.1 val- 
ues. We see here that the new calculation is now con- 
sistent with the low-density l imit . Meanwhile , the lower 
value found by lOsten et al.l (pOOl of i? = 2.79 ± 0.28 
(shown as the pink box) led to an inferred density at or 
near the low-density limit (less than 5 x lO^'^ cm~'^), while 
the new results show a measurable density from Ne IX in 
CrB, ranging from 5 x 10^° — 3.5 x 10^^ cm~^. 

Our results show that an update for the He-like isose- 
quence for all astrophysically important elements is ur- 
gently needed, since the existing data systematically un- 
derestimate key diagnostic line ratios. This process is 
underway, with new calculations being done by our and 
other groups. All available updates will be included in 
the next revision (v2.0.0) of the ATOMDB. 

We acknowledge support from the Chandra GO The- 
ory program. Support for NSB was provided by NASA 
contract NAS8-03060 to the Smithsonian Astrophysical 
Observatory for the CXC. 
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Fig. 3. — [Left] Ne IX R-ratio as a function of density at two different temperatures. [Rigfit] Ne IX R-ratio as a function of temperature 
at a density of 1 cm""^. 
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Fig. 4. — [Left] Ne IX G-ratio with horizontal lines showing the measurement of G from Capella. The best-fit temperature is logT = 
6.46 ± 0.14, in much better agreement with the temperature of maximum emissivity for Ne IX (log T = 6.6). [Right] Ne IX R-ratio plotted 
for all four calculations, at T = lO'^'^K and T = lO^ ^K. Shown in yellow is the range of R measured from Capella IjNess et al.ll2003bl ). 
while the pink box shows the measurement from CrB by lOsten et al.l l|2003l ) . 
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